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We calculate the divergences of the generating functional of quenched Chi- 

KJi \ ral Perturbation Theory to one loop for a generic number of flavours. The 

^ ■ flavour number dependence of our result enlightens the mechanism of quark 



loop cancellation in the quenched effective theory for any Green function or S 
matrix element. We also apply our results to tttt scattering and evaluate the 
coefficient of the chiral log in the S-wave scattering lengths for the quenched 
case. 
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Lattice simulations are at present the only available method to study in a quan- 
titative way the nonperturbative regime of Quantum Chromodynamics. Although 
theoretically well founded, the method is limited by a number of technical prob- 
lems, and by the computing power which is not yet enough to fully cope with QCD, 
despite continuous progress. One of the main technical difficulties is the implemen- 
tation of an efficient algorithm to evaluate the path integral over fermions, which 
is at present very time consuming. The usual way to go around this problem has 
been the use of the "quenched" approximation, i.e. setting to one the determinant 
resulting from the integration over fermions. In the language of Feynman diagrams, 
this approximation corresponds to neglecting quark loops. 

Quenching is so widely used that it has become important to understand a priori 
and in a quantitative way what is the effect introduced by this approximation |jl]]. A 
very promising approach to this problem is, in our opinion, the effective Lagrangian 
method proposed by Bernard and Golterman 0. This is based on the observation 
due to Morel (see Ref. |]l|) that the quenched approximation to QCD can be de- 
scribed by a Lagrangian where one adds to the usual quark fields a corresponding 
number of ghost (complex) scalar fields, with the same kinetic term, and the same 
coupling to gluons. The ghost fields exactly cancel the fermionic determinant. When 
the quark masses are set to zero the quenched QCD Lagrangian has a larger sym- 
metry than the usual U{N)l x U{N)ji, and must have a special fermionic character, 
since it transforms bosons into fermions. This extended symmetry is described by 
the graded group U{N\N)l x U{N\N)ji. Following exactly the same steps that lead 
from the globally symmetric QCD Lagrangian in the chiral limit to Chiral Pertur- 
bation Theory (CHPT) ||^, ^, Bernard and Golterman have proposed an effective 
Lagrangian that should describe the low energy regime of quenched QCD: quenched 
CHPT (often qCHPT in the following). The basic new ingredient is the presence 
of ghost fermion and boson fields that enter only inside loops. These are Goldstone 
particles of the extended symmetry, and can be viewed as bound states formed 
by a quark and a ghost-antiquark, or viceversa, and ghost-quark ghost-antiquark. 
Also for these new fields the structure of the interaction is dictated by the graded 
symmetry, which however leaves the coupling constants unconstrained. The main 
advantage of this method is that one has in principle a systematic way to study the 
modifications introduced by quenching: much as in CHPT, one can e.g. evaluate 
unambiguously the chiral logs present in any observable quantity. 

Quenched CHPT has been successful in showing that the chiral logs in the two 
point function of the axial current to one loop disappear (i.e. both in F^ and M^, 
see Ref. [^), and that the anomalous rj' field is ill-defined, developing a double pole 
in the propagator, as one would expect. The systematics of the cancellation of chiral 



logs for two point functions has been further investigated in Ref. |0. However, a 
general and strong argument in favour of the validity of this effective Lagrangian 
method is still missing. Such an argument could be produced if one were able to 
prove that the cancellation mechanism between pion and ghost loops is such that 
the pion loops that do not contain quark loops survive it (e.g. in vrvr scattering the 
graph in Fig. |l|a should remain whereas that in Fig. |lp3 not). This should be done 
in general, for any Green function or S matrix element. 

Although at first sight this may seem difficult, there is in fact a straightforward 
way to do it. The idea is to develop a path integral formulation of quenched CHPT 
much on the same line as that of ordinary CHPT. There, Gasser and Leutwyler Q] 
have shown that, by using the background field method and heat kernel techniques, 
it is possible to derive in closed form the divergent part of the one loop generating 
functional. This is a chiral invariant Lagrangian at order p^. The calculation was 
done in a general SU{N)ji x SU{N)l theory, so that the A^ dependence of the 
divergent part of the one loop generating functional is explicit. In the result of 
Gasser and Leutwyler four different powers of N occur: A^ with k = —2, —1,0, 1. 
Let us disregard for the moment the negative powers - we shall discuss them later. 
If we consider the case of A^ degenerate quarks, the interpretation of the terms 
linear in A^ and constant is rather straightforward: they come from the pion loops 
that contain one quark loop or none at all, respectively. This conclusion is based 
on the simple observation that since the strong interaction is flavour-blind, each 
quark loop produces a factor A^. Then, if quenched CHPT is working properly, the 
fermionic ghost loops have to cancel only the pion loops that give rise to a divergence 
proportional to A^. Our main aim is to verify that this is exactly what happens. 

For the sake of clarity we are going to ignore the presence of the U{1)a anomaly 
of QCD, in the following. Taking it into account would produce, as a main effect, a 
double pole in the r]' 0] . How to cope with the diseases of this part of the theory 
has been extensively analyzed by Bernard and Golterman 0, ^ looking at various 
different observables. In the perturbative treatment of the path integral we are 
considering, these contributions simply add to the standard pion loop contributions, 
and we can safely put them aside for the moment |^. 

The starting point of our calculation is the lowest order qCHPT Lagrangian, a 
graded symmetry generalization of the CHPT Lagrangian. This can be written as 



where str stands for the supertrace, the generalization of the trace to graded matrices 
(see Ref. H for details), D^^Us = d^^Us - ir^U^ + iUsl'^, Xs = 25o(s, + ip,), and 





the field Us can be written as Us = exp(V2i/F^), where $ is the hermitian non 
traceless block-matrix 

$ = 

containing the physical pseudoscalar field 0, the ghost field 0, both of bosonic nature, 
and the hybrid fields 9, 9^ of fermionic nature. The scalar field Ss contains the quark 
mass matrix Ai through: 

/ \/i n \ 

+ 6ss 

and that for our purposes is taken diagonal: A4 = niql. All the Goldstone bosons 
will then have the same mass: M^ = 2BQmq. The Lagrangian in Eq. (^ contains 
the external fields r^, l^, Ss, Ps, which are generalizations of the standard external 
fields, in order to make the Lagrangian locally invariant under the group U{N\N)ji x 
U{N\N)l. However, since we are not interested in Green functions of the spurious 
fields, we directly set to zero all the corresponding spurious external fields. The 
generating functional will then become a function of the usual external fields only: 
f^^i = Vfj, + a^, l^ = v^ — a^, s, p, which are all A^ x A^ matrices. 

To calculate quantum corrections we expand the leading order action given by 
the Lagrangian in Eq. (j^) in the vicinity of the classical solution which is de- 
termined by the external sources through the classical equations of motion. We 
define the classical solution as Us = u^ and describe the fiuctuations around it as 
Us = Usexp{i^s)us = ^s(l + iC,s — 1/2^^ + • • •)■"«• To get the generating functional 
to one loop we need to expand the action in the ^s field up to second order. Skip- 
ping all the details [^ we end up with the following gaussian integral (in Minkowski 
space-time): 



detD 
{detDf 



•^ .. ,^,1/2 > (2) 



where we used the decomposition of the fields ^,^, C)C^ (quantum fiuctuations cor- 
responding to the classical fields (p, (p, 9, 9\ respectively) in terms of the A^^ 
generators of each U{N) fiavour subgroup as ,^ = C'Xa, with \a = Aa/v2 for 
a = 1, . . . , A^^ — 1 (the A^'s are the usual Gell-Mann matrices of SU{N)), and 
Aq = 1/vN. We remark that the scalar ghosts ^ decouple from the physical pions 
and the integral over them produces only an irrelevant constant. We also stress 



that both differential operators D"^'' and D"'' are now functions only of the standard 
external fields and the field at the classical solution, since we have put to zero all 
the spurious external fields. The differential operator D"^ is defined as follows 



d,C = d,C + Tt^b , 



where 



'^ab 



-(r,[r,A^]) , a'^'' = ^(K,r]K,A^])-i({r,A^}x+) . 



(3) 



(4) 



Let us recall some standard definitions in CHPT: F^ = l/2{[u\dfj_u\ — iu^r^u — 
iul^u^), Ufj_ = iv)D^Uu\ and x+ = w^xm^ + ux'^u. The differential operator D"'^ 
acting on the ghost field ( is defined like in Eq. (H), but with barred quantities, 
given by 

1 



(r^A'^A") 



-ab 



^{{uX + X+ + ^BoM)X^X'') , 



(5) 



where Ad is again the quark mass matrix, which is also contained in the external 
scalar field s = A4 + 6s. 

The divergent part of Eq. (0) can be derived in closed form by regularizing the 
determinants in d dimensions and using standard heat kernel techniques. The result 
reads: 

_±_/,,|^(r„„r-) + l 



In det D 



l(«,n.)K«-) + ^(nX)' 



+ y ((wX)') + j(wXx+) + -(mX)(x+) 



+ 



i In det D 



-1 



(47r)2(rf-4) 
and their difference gives 
-1 



dx 



N N 

D Id 



(6) 

+ ••(7) 



r^qCHPT 

one loop 



(1 1 1 

-^M'{uX)-jM\x+)} + ... . (8) 



The ellipses stand everywhere for the finite contributions to the one loop generating 
functional. Eq. (^) shows explicitly the fiavour dependence of the qCHPT functional 
to one loop and the comparison with Eq. (^ provides the result we were after. We 
stress a few important points: 

1. Last line of Eq. (^) contains extra contributions with respect to the SU{N) x 
SU{N) determinant of Gasser and Leutwyler |Q, due to the presence of the 
singlet field. Two of those terms are proportional to {u^) = V^0o; and only 
contribute to processes with external singlet fields. The last two terms arise 
from loops with the singlet field running inside. These terms exactly remove 
the negative power dependence upon N in the rest of the expression. This 
shows that the presence of such terms in the SU{N) theory can be understood 
as an effect of the U{1)a anomaly that decouples the singlet field. 

2. The ( field loop in Eq. (^ produces only contributions linear in A^ as ex- 
pected, and cancels completely the terms linear in A^ of the full generating 
functional. We also note that the N dependence in Eq. (|^) is not fully ex- 
plicit. In fact, since the expansion of x+ in powers of starts with a constant 
term proportional to the quark mass matrix, its trace is proportional to A^ in 
the degenerate case we are considering here. The last two terms in Eq. (H) 
cancel this dependence in the final expression. This result shows that quenched 
CHPT does what it should do, i.e. it contains only the pion loops that do not 
contain quark loops. 

3. Each pole at c? = 4 produces a chiral logarithm after renormalization. A 
different way to present our result is to say that we have calculated the chiral 
logs of any Green function in qCHPT to one loop. We obtained that these 
are nonzero in general, though certainly different from those of the ordinary 
CHPT case. How different has to be investigated case by case, since it is not 
possible to identify a general behaviour from our result. 

In particular, from Eq. (^) it is easy to verify the results already obtained by Bernard 
and Golterman about the absence of chiral logs in the chiral condensates, pion 
masses and decay constants. The situation is different when one considers Green 
functions with more than two external legs. To see what happens in one concrete 
example, we consider the tttt scattering amplitude - also because there are lattice 
calculations of the two S'-wave scattering lengths available, Ref. 0. As we argued 
above, the presence of chiral logs even in the quenched theory has to be interpreted 
as due to diagrams with pion loops that do not contain quark loops. For the vrvr 
scattering amplitude an example is given in Fig. |l[ 





(a) 



(b) 



Figure 1: Two examples of pion loop graphs contributing to tttt scattering in the 
quark-flow diagram picture (all lines are quark lines). Diagram (a) does not contain 
quark loops, whereas diagram (b) does. 



The complete amplitude in quenched CHPT reads: 

s-Ml 1 r 1 / . 1 



Afr^is^t^u) 



4'^ 



Ss^ + (t - uf 

2 -, 



F2 4F4 I 2 V 167r2 

+ s2J(s) + (t - 2M2)' Jit) + {u- 2MlY J{u 
+ c^Mt + C2sMl + cl{^l)s^ + cl{^l){t - uf) 
+ 5'(^o,a)(s,t,M) + 0(/) , 



where J{q^) = J{q^) - J(0), with 



J{q' 



dH 



1 



I J {2'nY (M2 - P) (M2 - (/ - g)2) 



(9) 



(10) 



and L = (IGtt^) ^ log(M^//i^j. The polynomial part with coefficients ci,...C4 
comes from the counterterm Lagrangian at order p^. «S'(mo,o)(s, t, m) is the renormal- 
ized contribution to the amplitude of the singlet loops with one {itiq and a) vertex 
insertion on the singlet propagator. The effects and diseases of S(mo,a){.s,t,u) have 
been discussed in |^, |TD|. F^, and M^ are the renormalized quenched values of F 
and M at order p^. While F remains unchanged to one loop, M gets a one loop 
correction from {itlq and a) vertex insertion which we included in S'(mo,a)('S, t,M). 
Finite contributions to both F^ and M^^ from the order p^ Lagrangian are meant to 
be included in the renormalized polynomial part. Let us stress that the divergent 
part of the amplitude before renormalization, and the corresponding chiral logs in 
the renormalized amplitude Eq. (P) can also be obtained from the generating func- 
tional, Eq. (^), by expanding it in powers of the external fields for N = 2, and 
taking the coefficient of the relevant term. This offers a welcome check on the full 
calculation. 



We focus now on the S-wave scattering lengths. In full CHPT it is well known 



that the chiral logs dominate the one loop correction at /x = 1 GeV |11] . In the 
quenched case, the coefficients of the chiral logs are: 

«S^^"^^ = ^.i^ + ^(-fL + ..)+OiM-)\ , (11) 




2 



^^^"^^ = -J^2[^ + ^(^L + ...) + 0{Mt)\ , (12) 

while in full CHPT the corresponding coefficients are: —9/2 for / = and 3/2 for 




1 = 2 ||TT|. In both cases the change in the correction due to quenching is of the order 
of 30%; however the one loop correction itself is roughly of this size in full CHPT (see 
Ref. |[n[]), so that the overall relative change can be estimated to be around 10%. 
Of course this is not the whole story: we have not included the remaining analytic 
contributions from the regular part of the amplitude, the constants Cj's in Eq. (^, 
and the contributions from the ill part of the amplitude, S(^rno,a){s,t,u). As for the 
first point, if one assumes that quenching will not change the order of magnitude 
of the finite part of the 0{p'^) constants at /i = 1 GeV (one could argue that some 
kind of Vector Meson Dominance should be valid also in the quenched case), one 
can conclude that those contributions should still produce a small correction. 

The ill part of the amplitude is in principle more dangerous and requires a more 
careful treatment, as was done by Bernard and Golterman 0. They have shown in 
fact that this part dramatically changes the very method to compute the scattering 
lengths on the lattice: the Liischer's formula |T^ that relates finite volume effects 
to the S-wave scattering lengths is modified by quenching. However, the numerical 
analysis of this effect, adapted to the staggered fermions calculation by Fukugita 
et al. IP (the only calculation done with a pion mass small enough to justify the 
application of CHPT) indicates that this effect should be small. Putting the various 
pieces together, we can conclude that a complete analysis in qCHPT to one loop of 
the TTTT scattering amphtude, seems to provide an explanation for the rather good 
agreement between the lattice calculation of S-wave scattering lengths done in Ref. 
and standard CHPT, as was pointed out by one of us in Ref. [^ . 



Let us summarize our results: we have calculated the divergent part of the 
generating functional of quenched CHPT to one loop, Eq. (|]), and found that, for 
degenerate quark masses, it does not contain any explicit flavour number dependence. 
Arguing that for a non-anomalous theory this dependence may arise only through 
the presence of quark loops (become manifest here through pion loops), we conclude 
that quenched CHPT is working properly. In our opinion this result puts on a 
very solid basis quenched CHPT, the effective theory proposed by Bernard and 
Golterman M to describe the low energy physics of quenched QCD. Our calculation 



provides also the chiral logs of any Green function, and shows that these are nonzero 
in general. As an example we have discussed the case of the vrvr scattering amplitude 
and found that the coefficient of the chiral log of the two S'-wave scattering lengths 
has been modified by 30%, a rather modest effect. 
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